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Abstract. On 23"^"^ November 2009, during the early commissioning of the CERN Large Hadron Collider 
(LHC), two counter- rotating proton bunches were circulated for the first time concurrently in the machine, 
at the LHC injection energy of 450 GeV per beam. Although the proton intensity was very low, with 
only one pilot bunch per beam, and no systematic attempt was made to optimize the collision optics, all 
LHC experiments reported a number of collision candidates. In the ALICE experiment, the collision region 
was centred very well in both the longitudinal and transverse directions and 284 events were recorded in 
coincidence with the two passing proton bunches. The events were immediately reconstructed and analyzed 
both online and offline. We have used these events to measure the pseudorapidity density of charged primary 
particles in the central region. In the range |77| < 0.5, we obtain diVch/dr? = 3.10±0.13(staf.) ±0.22{syst.) 
for all inelastic interactions, and dA^ch/d»? = 3.51 ± 0.15{stat.) ± 0.25{syst.) for non-single diffractive 
interactions. These results arc consistent with previous measurements in proton-antiproton interactions at 
the same centre-of-mass energy at the CERN SppS collider. They also illustrate the excellent functioning 
and rapid progress of the LHC accelerator, and of both the hardware and software of the ALICE experiment, 
in this early start-up phase. 
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1 Introduction 

The very first proton-proton collisions at Point 2 of the 
CERN Large Hadron Collider (LHC) [l] occurred in the 
afternoon of 23'''^ November 2009, at a centre-of-mass en- 
ergy = 900 GeV, during the commissioning of the ac- 
celerator. This publication, based on 284 events recorded 
in the ALICE detector on that day, describes a deter- 
mination of the pseudorapidity density of charged primary 
particlesj^ diVch /dry (r/ = — In tan 0/2, where 9 is the polar 
angle with respect to the beam line) in the central pseudo- 
rapidity region. The purpose of this study is to compare 
with previous measurements for proton-antiproton (pp) 
collisions at the same energy [3^, and to establish a ref- 
erence for comparison with forthcoming measurements at 
higher LHC energies. 

The event sample collected with our trigger contains 
three different classes of inelastic interactions, i.e. colli- 
sions where new particles are produced: non-diffractive, 
single-diffractive, and double-diffractiv^ Experimentally 

^ Here, primary particles are defined as prompt particles pro- 
duced in the collision and all decay products, except products 
from weak decays of strange particles such as and A. 

^ Inelastic pp collisions are usually divided into these classes 
depending on the fate of the interacting protons. If one (both) 



we cannot distinguish between these classes, which, how- 
ever, are selected by our trigger with different efficiencie^ 

In order to compare our data with those of other exper- 
iments, we provide the result with two different normal- 
izations: the first one (INEL) corresponds to the sum of all 
inelastic interactions and corrects the trigger bias individ- 
ually for all event classes, by weighting them, each with its 
own estimated trigger efficiency and abundance. The sec- 
ond normalization (non-single-diffractive or NSD) applies 
this correction for non-diffractive and double-diffractive 
processes only, while removing, on average, the single- 
diffractive contribution. 

Multiparticle production is rather successfully describ- 
ed by phenomenological models with Pomeron exchange, 
which dominates at high energies pllsl. These models re- 



incoming beam particle(s) are excited into a high-mass state, 
the process is called single (double) diffraction; otherwise the 
events are classified as non-diffractive. Particles emitted in 
diffractive reactions are usually found at rapidities close to 
that of the parent proton. 

We estimate the trigger efficiency for each class using the 
process-type information provided by Monte Carlo generators; 
the values vary by up to a factor of two between classes and 
are listed in Section 3. The relative abundance of each class is 
taken from published data (see text). 
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late the energy dependence of the total cross section to 
that of the multiplicity production using a small num- 
ber of parameters, and are the basis for several Monte 
Carlo event generators describing soft hadron collisions 
(see for example [6||8]). According to these models, it is 
expected that the charged-particle density increases by a 
factor 1.7 and 1.9 when raising the LHC centre-of-mass 
energy from 900 GeV to 7 and 14TeV respectively (i.e. 
intermediate and nominal LHC energies). The difference 
in charged-particle densities between pp and pp interac- 
tions is predicted to decrease as 1/ at high energies j9| . 
This difference was last measured at the CERN ISR to be 
in the range 1.5-3% flO] at ^/s = 53 GeV. Extrapolating 
these values to ^ys — 900GeV, one obtains a very small 
difference of about 0.1-0.2%. Therefore, we will compare 
our measurement to existing pp data and also to different 
Monte Carlo models. 

This article is organized as follows: Section 2 describes 
the experimental conditions during data taking; the main 
features of the ALICE detector subsystems used for this 
analysis are decribed in Section 3; Section 4 is dedicated 
to the event selection and data analysis; the results are 
discussed in Section 5 and Section 6 contains the conclu- 
sion. 



2 LHC and the run conditions 

The LHC, built at CERN in the circular tunnel of 27 km 
circumference previously used by the Large Electron-Posi- 
tron collider (LEP), will provide the highest energy ever 
explored with particle accelerators. It is designed to col- 
lide two counter-rotating beams of protons or heavy ions. 
The nominal centre-of-mass energy for proton-proton col- 
lisions is 14TeV. However, collisions can be obtained down 
to ^/s = 900 GeV, which corresponds to the beam injec- 
tion energy. 

The results from the first proton-proton collisions pre- 
sented here were obtained during the early commissioning 
phase of the LHC, when two proton bunches were circu- 
lating for the first time concurrently in the machine. The 
bunches used were the so-called "pilot bunches" : low in- 
tensity bunches used during machine commissioning, with 
a few 10^ protons per bunch. The two beams were brought 
into nominal position for collisions without a specific at- 
tempt to maximize the interaction rate. The nominal r.m.s. 
size of LHC beams at injection energy is about 300 /im in 
the transverse direction and 10.5 cm in the longitudinal 
(z-axis) direction. However, at this early stage, the beam 
parameters can deviate from these nominal values; they 
were not measured for the fill used in this analysis. For 
the previous fill, for which the longitudinal size was mea- 
sured, it was found to be shorter, with an r.m.s. of about 
8 cm. Assuming Gaussian beam profiles, the luminous re- 
gion should be smaller than the beam size by a factor of 
a/2 in all directions. 

Shortly after circulating beams were established, the 
ALICE data aquisition system [TT] started collecting ev- 
ents with a trigger based on the Silicon Pixel Detector 



(SPD), requiring two or more hits in the SPD in coinci- 
dence with the passage of the two colliding bunches as in- 
ferred from beam pickup detectors. As a precaution, only 
a small subset of the detector subsystems, including the 
silicon tracking detectors and the scintillator trigger coun- 
ters, was turned on, in order to assess the beam conditions 
provided by the LHC. 

The trigger rate was measured just before collisions 
with the same trigger conditions. Without beams we mea- 
sured a rate of 3 x 10"'* Hz (in coincidence with one bunch 
crossing interval per orbit). In coincidence with the pas- 
sage of the bunch of one circulating beam the rate was 
0.006 Hz. As soon as the second beam was injected in 
the accelerator, the event rate increased significantly, to 
0.11 Hz. The first event that was analyzed and displayed 
in the counting room by the offline reconstruction software 
AliRoot [12] running in online mode is shown in Fig. [l] 
This marked symbolically the keenly anticipated start of 
the physics exploitation of the ALICE experimeni]^ The 
online reconstruction software implemented in the High- 
Level Trigger (HLT) computer farm [l3] also analyzed the 
events in real time and calculated the vertex position of 
the collected events, shown in Fig. 2j The distributions 
are very narrow in the transverse plane (sub-millimetre, 
including contributions from detector resolution and resid- 
ual misalignment), of about the expected size in the lon- 
gitudinal direction and well positioned with respect to the 
nominal centre of the ALICE detector. This provided im- 
mediate evidence that a substantial fraction of the events 
corresponded to collisions between the protons of the two 
counter-rotating beams. 

After 43 minutes, the two beams were dumped in order 
to proceed with the LHC commissioning programme. In 
total, 284 events were triggered and recorded during this 
short, but important, first run of the ALICE experiment 
with colliding beams. 



3 The ALICE experiment 

ALICE, designed as the dedicated heavy-ion experiment 
at the LHC, also has excellent performance for proton- 
proton interactions [14^ . The experiment consists of a large 
number of detector subsystems [2] inside a solenoidal mag- 
net {B = 0.5 T). The magnet was off during this run. 

During the several months of running with cosmic rays 
in 2008 and 2009, all of the ALICE detector subsystems 
were extensively commissioned, calibrated and used for 
Data were collected for an initial align- 



data taking 15 



ment of the parts of the detector that had sufficient ex- 
posure to the mostly vertical cosmic ray flux. Data were 
also taken during various LHC injection tests to perform 
timing measurements and other calibrations. 

Collisions take place at the centre of the ALICE de- 
tector, inside a beryllium vacuum beam pipe (3 cm in ra- 
dius and 800 fim thick) . The tracking system in the AL- 
ICE central barrel covers the full azimuthal range in the 



* The event display started shortly after data taking and 
therefore missed the first few events. 
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Fig. 1. The first pp collision candidate shown by the event display in the ALICE counting room (3D view, r-ij) and r-z 
projections), the dimensions are shown in cm. The dots correspond to hits in the silicon vertex detectors (SPD, SDD and SSD), 
the lines correspond to tracks reconstructed using loose quality cuts. The ellipse drawn in the middle of the detector surrounds 
the reconstructed event vertex. 
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Fig. 2. Online display of the vertex positions reconstructed by the High-Level Trigger (HLT). The figure shows, counter- 
clockwise from top left, the position in the transverse plane for all events with a reconstructed vertex, the projections along the 
transverse coordinates x and and the distribution along the beam line (2-axis). 
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Fig. 3. Arrival time of particles in the VZERO detectors relative to the beam crossing time (time zero). A number of beam- 
halo or beam-gas events are visible as secondary peaks in VZERO- A (left panel) and VZERO-C (right panel). This is because 
particles produced in background interactions arrive at earlier times in one or the other of the two counters. The majority of 
the signals have the correct arrival time expected for collisions around the nominal vertex. 



pseudorapidity window \r]\ < 0.9. It has been designed to 
cope with the highest charged-particle densities expected 
in central Pb-Pb collisions. The following four detector 
subsystems were active during data taking and were used 
in this analysis. 

— The Silicon Pixel Detector (SPD) consists of two cylin- 
drical layers with radii of 3.9 and 7.6 cm and has about 
9.8 million pixels of size 50 x 425 /im^. It covers the 
pseudorapidity ranges \ri\ < 2 and I77I < 1.4 for the 
inner and outer layers respectively, for particles orig- 
inating at the centre of the detector. The effective rj- 
acceptance is larger due to the longitudinal spread of 
the position of the interaction vertex. The detector is 
read out by custom-designed ASICs bump-bonded di- 
rectly on silicon ladders. Each chip contains 8192 chan- 
nels and also provides a fast trigger signal if at least 
one of its pixels is hit. The trigger signals from all 
1200 chips are then combined in a programmable logic 
unit which provides a level-0 trigger signal to the cen- 
tral trigger processor. The total thickness of the SPD 
amounts to about 2.3% of a radiation length. About 
83 % of the channels were operational for particle de- 
tection and 77 % of the chips were used in the trigger 
logic. The SPD was aligned using cosmic-ray tracks 
collected during 2008 [16], and the residual misalign- 
ment was estimated to be below 10 /im for the modules 
well covered by mostly vertical tracks. The modules 
on the sides are likely to be affected by larger residual 
misalignment. 

— The Silicon Drift Detector (SDD) consists of two cylin- 
drical layers at radii of 15.0 and 23.9 cm and covers 
the region \ri\ < 0.9. It is composed of 260 sensors 
with an internal voltage divider providing a drift field 
of 500 V/cm and MOS charge injectors that allow 
measurement of the drift speed via dedicated calibra- 
tion triggers. The charge signal of each of the 133 000 
collection anodes, arranged with a pitch of 294 fiia, 



is sampled every 50 ns by an ADC in the front-end 
electronics. The total thickness of the SDD layers (in- 
cluding mechanical supports and front-end electronics) 
amounts to 2.4 % of a radiation length. About 92 % of 
the anodes were fully operational. 

— The two layers of the double-sided Silicon Strip Detec- 
tor (SSD) are located at radii of 38 and 43 cm respec- 
tively, covering I??] < 0.97. The SSD consists of 1698 
sensors with a strip pitch of 95 fim and a stereo angle 
of 35 mrad. The detector provides a measurement of 
the charge deposited in each of its 2.5 x 10^ strips. The 
position resolution is better than 20 /im in the r-ip di- 
rection and about 0.8 mm in the direction along the 
beam line. The thickness of the SSD, including sup- 
ports and services, corresponds to 2.2 % of a radiation 
length. About 90% of the SSD area was active during 
data taking. 

— The VZERO detector consists of two arrays of 32 scin- 
tillators each, which are placed around the beam pipe 
on either side of the interaction region: VZERO- A at 
z = 3.3 m, covering the pseudorapidity range 2.8 < 
1] < 5.1, and VZERO-C at z = —0.9 m, covering the 
pseudorapidity range —3.7 < j] < —1.7. The time reso- 
lution of this detector is better than 1 ns. Its response 
is recorded in a time window of ±25 nsec around the 
nominal beam crossing time. For events collected in 
this run, the arrival times of particles at the detector 
relative to this "time zero" is shown in Fig. [3] Note 
that in general several particles are registered for each 
event. Particles hitting one of the detectors before the 
beam crossing have negative arrival times and are typ- 
ically due to interactions taking place outside the cen- 
tral region of ALICE. 

More details about the ALICE experiment and its detector 
subsystems can be found in 

The trigger used to record the events for the present 
analysis is defined by requiring at least two hit chips in 
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the SPD, in coincidence with the signals from the two 
beam pick-up counters indicating the presence of two pass- 
ing proton bunches. The efficiency of this trigger as well 
as all other corrections have been studied using two dif- 
ferent Monte Carlo generators, PYTHIA 6.4.14 [l7] tune 
D6T [18] and PHOJET [s], for INEL and NSD interac- 
tions. The trigger efficiencies for non-diffractive, single- 
diffractive, and double-diffractive events were evaluated 
separately, and found to be 98-99%, 48-58%, and 53- 
76 % respectively. The ranges are determined by the two 
event generators. These event classes were combined for 
the corrections using the fractions measured by UA5 |19j : 
non-diffractive 0.767 ± 0.059; single-diffractive 0.153± 
±0.031; double-diffractive 0.08 ± 0.05. The resulting ef- 
ficiencies were found to be 87-91 % for the INEL normal- 
ization and 94-97% for the NSD normalization, again de- 
pending on the event generator used. 

The results presented in the following sections are those 
obtained with PYTHIA. The difference between results 
corrected with PYTHIA and PHOJET is used in the es- 
timate of the systematic uncertainty. 
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Fig. 4. Longitudinal vertex distribution from hit correlations 
in the two pixel layers of the ALICE inner tracking system. 
Vertical dashed lines indicate the region \z\ < 10 cm, where 
the events for the present analysis are selected. A Gaussian fit 
with an estimated r.m.s. of about 4 cm to the central part is 
also shown. 



4 Data analysis 

The data sample used in the present analysis consists of 
284 events recorded without magnetic field. The results 
presented here are based on the analysis of the SPD data. 
However, information from the SDD, SSD and VZERO 
was used to crosscheck the identification and removal of 
background events. 

In the SPD analysis, the position of the interaction 
vertex is reconstructed 
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by correlating hits in the two 
silicon-pixel layers to obtain tracklets. The achieved reso- 
lution depends on the track multiplicity and for this spe- 
cific vertex reconstruction is approximately 0.1-0.3 mm in 
the longitudinal direction and 0.2-0.5 mm in the trans- 
verse direction. For events with only one charged track, 
the vertex position is determined by intersecting the SPD 
tracklet with the mean beam axis determined from the 
vertex positions of other events in the sample. A vertex 
was reconstructed in 94% of the selected events. The dis- 
tribution of the vertex position in the longitudinal direc- 
tion (z-axis) is shown in Fig. W| For events originating 
from the centre of the detector, the vertex-reconstruction 
efficiency was estimated, using Monte Carlo simulations, 
to be 84 % for INEL interactions and 92 % for NSD col- 
lisions. These efficiencies decrease for larger jzj-values of 
the vertex in low-multiplicity events; therefore, only events 
with vertices within \z\ < 10cm were used. This allows for 
an accurate charged-particle density measurement in the 
pseudorapidity range \t]\ < 1.6 using both SPD layers. 

Using the reconstructed vertex as the origin, we calcu- 
late the differences in azimuthal {Aip, bending plane) and 
polar {A6, non-bending direction) angles of pairs of hits 
with one hit in each SPD layer. These tracklets [21] are 
selected by a cut on the sum of the squares of A(p and A9, 
each normalized to its estimated resolution (80 mrad and 
25 mrad, respectively) . When more than one hit in a layer 
matches a hit in the other layer, only the hit combination 



with the smallest angular difference is used. This occurs 
in only 2 % of the matched hits. 

The number of primary charged particles is estimated 
by counting the number of tracklets. This number was 
corrected for: 

— trigger inefficiency; 

— detector and reconstruction inefficiencies; 

— contamination by decay products of long-lived parti- 
cles (Kg, A, etc.), gamma conversions and secondary 
interactions. 

The corrections are determined as a function of the z- 
position of the primary vertex, and on the pseudorapidity 
of the tracklet. For the analyzed sample the average cor- 
rection factor for tracklets is about 1.5. 

The beam-gas and beam-halo background events were 
removed by a cut on the ratio between the number of 
tracklets and the total number of hits in the tracking 
system (SPD, SDD, and SSD); this ratio is smaller for 
background events (as measured in the previous fills trig- 
gering on the bunch passage from one side) than for col- 
lisions 
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In addition, the timing information from the 
VZERCTaetector was used for background rejection by re- 
moving events with negative arrival time (see Fig. pi. The 
event quality and event classification was crosschecked by 
a visual scan of the whole event sample. In total 29 events 
(i.e. about 10%) were rejected as beam induced back- 
ground, which is consistent with the rate expected from 
previous fills. The remaining background was estimated 
from the vertex distribution and found to be negligible. 
The contamination from coincidence with a cosmic event 
was estimated to be one event in the full sample. Indeed, 
two cosmic events were identified by scanning, both with- 
out reconstructed vertex. 

Particular attention has been paid to events having 
zero or one charged tracklets in the SPD acceptance. The 
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Table 1. Contributions to systematic uncertainties on the 
measurement of the charged-particle pseudorapidity density. 



2 4 6 8 10 

Generated multiplicity in |t|| < 1.4 

Fig. 5. Multiphcity dependence of the combined efficiency 
to select an event as minimum bias and to reconstruct its 
vertex in SPD, for non-diffractive (crosses), single-difi'ractive 
(squares), and double-diffractive (circles) events, based on 
PYTHIA events. 



vertex-finding efficiency for events with one charged par- 
ticle in the acceptance is about 80 %. The number of zero- 
track events has been estimated by Monte Carlo calcula- 
tions. The total number of collisions used for the normal- 
ization was calculated from the number of events selected 
for the analysis, corrected for the vertex-reconstruction in- 
efficiency. In order to obtain the normalization for INEL 
and NSD events, we further corrected the number of se- 
lected events for the trigger efficiency for these two event 
classes. In addition, for NSD events, we subtract the single- 
diffractive contribution. These corrections, as well as those 
for the vertex finding efficiency, depend on the event char- 
ged-particle multiplicity, see Fig. |5] The dependence of 
the event-finding efficiency (combining event selection and 
vertex finding) on multiplicity was calculated for differ- 
ent interaction types using our detector simulation, and is 
above 98% for events with at least two charged particles. 
The averaged combined corrections for the vertex recon- 
struction efficiency and the selection efficiency is 20 % for 
INEL interactions and much smaller for NSD interactions, 
due to the cancelation of some contributions. 

The various corrections mentioned above were calcu- 
lated using the full GEANT 3 i23] simulation of the AL- 
ICE detector as included in the offline framework Ali- 
Root. In order to estimate the systematic uncertainties, 
the above analysis was repeated by: 



Uncertainty 


Tracklet selection cuts 


need. 


Material budget 


negl. 


Misalignment 


0.5% 


Particle composition 


negl. 


Transverse-momentum spectrum 


0.5% 


Contribution of difi'raction (INEL) 


4% 


Contribution of difi'raction (NSD) 


4.5% 


Event-generator dependence (INEL) 


4% 


Event-generator dependence (NSD) 


3% 


Detector efficiency 


4% 


SPD triggering efficiency 


2% 


Background events 


negl. 


Total (INEL) 


7.2% 


Total (NSD) 


7.1% 



— evaluating the uncertainty in the normalization to 
INEL and NSD samples by varying the ratios of the 
non-diffractive, single-diffractive and double-diffractive 
cross sections according to their measured values and 
and using two different models for diffrac- 



errors 
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tion kinematics (PYTHIA and PHOJET). 

An additional source of systematic error comes from 
the limited statistics used so far to determine the efficien- 
cies of the SPD detector modules. In test beams, the SPD 
efficiency in active areas was measured to be higher than 
99.8%. This was crosschecked in-situ with cosmic data, 
but only over a limited area and with limited statistics. 
At this stage, we have assigned a conservative value of 4 % 
to this uncertainty. The triggering efficiency of the SPD 
was estimated from the data itself, using the trigger infor- 
mation recorded in the data stream for events with more 
than one tracklet, and found to be very close to 100%, 
with an error of about 2% (due to the limited statistics). 

These contributions to the systematic uncertainty on 
the charged particle pseudorapidity density are summa- 
rized in Table [l] Our conclusion is that the total system- 
atic uncertainty on the pseudorapidity density is less than 
±7.2 % for INEL coUisions and ±7.1 % for NSD coUisions. 
The largest contribution comes from uncertainties in cross 
sections of diffractive processes and their kinematic simu- 
lation. 

More details about this analysis, corrections, and the 
evaluation of the systematic uncertainties can be found 
in [24 . 



— applying different cuts for the tracklet definition (vary- 
ing the angle cut-off by ±50 %); 

— varying by ±10% the density of the material in the 
tracking system, thus changing the material budget; 

— using the non-aligned geometry; 

— varying by ±30 % the composition of the produced par- 
ticle types with respect to the yields suggested by the 
event generators; 

— varying the particle yield below lOOMeV/c by ±30%; 



5 Results 

Figure [6] shows the charged primary particle pseudorapid- 
ity density distributions obtained for INEL and NSD inter- 
actions in the range \r/\ < 1.6. The pseudorapidity density 
obtained in the central region \ri\ < 0.5 for INEL interac- 
tions is 3.10±0.13{stat.)±0.22{syst.) and for NSD interac- 
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Table 2. Comparison of charged primary particle pseudorapidity densities at central pseudorapidity (|7;| < 0.5) for inelastic 
(INEL) and non-single diffractive (NSD) collisions measured by the ALICE detector in pp interactions and by UA5 in pp 
interactions 3 at a centre-of-mass energy of 900 GeV. For ALICE, the first error is statistical and the second is systematic; no 
systematic error is quoted by UA5. The experimental data are also compared to the predictions for pp collisions from different 
models. For PYTHIA the tune versions are given in parentheses. The correspondence is as follows: D6T is tune (109); ATLAS 
CSC is tune (306); Perugia-0 is tune (320). 



Experiment 
Model 


ALICE pp 


UA5 pp D 


QGSM 


25 


(109) 


18 


PYTHI7 
(306) 


2^ 


] 


L7 


(320) 


27 


PHOJET [s] 


INEL 
NSD 


3.10 ±0.13 ±0.22 
3.51 ±0.15 ±0.25 


3.09 ±0.05 
3.43 ±0.05 


2.98 
3.47 




2.33 
2.83 


2.99 
3.68 




2.46 
3.02 


3.14 
3.61 




1- 



i 



n/s = 900 GeV 
■ ALICE pp NSD i UA5 pp NSD 
□ ALICE pp INEL A UASpplNEL 



-1 



1 



Fig. 6. Pseudorapidity dependence of dA^'ch/d?) for INEL and 
NSD collisions. The ALICE measurements (squares) are com- 
pared to UA5 data (triangles) [sj. The errors shown are statis- 
tical only. 



iSR (pp) INEL 
UA1 (pp) NSD 
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Fig. 7. Charged-particle pseudorapidity density in the central 
rapidity region in proton-proton and proton-antiproton inter- 
actions as a function of the centre-of-mass energy. The dashed 
and solid lines (for INEL and NSD interactions respectively) 
indicate the fit using a power-law dependence on energy. 



tions is 3.51±0.15(stot)±0.25(s?/st). Also shown in Fig.|6] 
are the previous measurements of proton-antiproton inter- 
actions from the UA5 experiment ^ . Our results obtained 
for proton-proton interactions are consistent with those 
for proton-antiproton interactions, as expected from the 
fact that the predicted difference (0.1-0.2 %) is well below 
measurement uncertainties. The measurements at central 
pseudorapidity {\rj\ < 0.5) are summarized in Table |2] 
together with model predictions obtained with QGSM, 
PHOJET and three different PYTHIA tunes. PYTHIA 
6.4.14, tune D6T, and PHOJET yield respectively the low- 
est and highest charged particle densities. Therefore, these 
two have been used for the evaluation of our systematic 
errors. PYTHIA 6.4.20, tunes ATLAS CSC and Perugia-0, 
are candidates for use by the LHC experiments at higher 
LHC energies and are shown for comparison. 

Figure |7] shows the centre-of-mass energy dependence 
of the pseudorapidity density in the central region (|?7| < 
0.5). The data points are obtained in th e [77I < .5 range 
from this experiment and from references 3|10|28 - 3l] , and 
are corrected for differences in pseudorapidity range where 
necessary, fitting the pseudorapidity distribution around 
77 = 0. As noted above, there is good agreement between 
pp and pp data at the same energy. The dashed and solid 



lines (for INEL and NSD interactions respectively) are 
obtained by fitting the density of charged particles in the 
central pseudorapidity rapidity region with a power-law 
dependence on energy. 

Using this parametrization, the extrapolation to the 
nominal LHC energy of -y/s — 14TeV yields dNch/di] = 
5.5 and diVch/dry = 5.9 for INEL and NSD interactions 
respectively. 



6 Conclusion 

Proton-proton collisions observed with the ALICE detec- 
tor in the early phase of the LHC commissioning have been 
used to measure the pseudorapidity density of charged pri- 
mary particles at y/s — 900 GeV. In the central pseudo- 
rapidity region {\ri\ < 0.5), we obtain dNch/drj = 3.10 ± 
0.13{stat.) ± 0.22{syst.) for all inelastic and dNch/drj — 
3.51 ± 0.15{stat.) ± 0.25{syst.) for non-single diffractive 
proton-proton interactions. The results are consistent with 
earlier measurements of primary charged-particle produc- 
tion in proton-antiproton interactions at the same energy. 
They are also compared with model calculations. 
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These results have been obtained with a small sam- 
ple of events during the early commissioning of the LHC. 
They demonstrate that the LHC and its experiments have 
finally entered the phase of physics exploitation, within 
days of starting up the accelerator complex in November 
2009. 
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